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ABSTRACT 


In the study of normal evaporation, it is assumed that the 
evaporating alloy is homogeneous, that the vapor is instantly 
removed, and that the alloy follows Raoult's law. The differen- 
tial equation of normal evaporation relating the evaporating time 
to the final solute concentration is given and solved for several 
important special cases . Uses of the derived equations are ex- 
emplified with a Ni-Al alloy and some binary iron alloys. The 
accuracy of the predicted results are checked by analyses of 
actual experimental data on Fe~Ni and Ni~Cr alloys evaporated 
at 1600 °C, and also on the vacuum purification of beryllium. 
These analyses suggest that the normal evaporation equations pre- 
sented here give satisfactory results that are accurate to within 
an order of magnitude of the correct values, even for some highly 
concentrated solutions. Limited diffusion and the resultant sur- 
face solute depletion or enrichment appear important in the ex- 
tension of this normal evaporation approach. 
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INTRODUCTION 


In normal evaporation, we assume that the evaporating alloy 
is always homogeneous in composition, that the vapor is instantly 
removed, and that the alloy follows the Raoul t’s law (Ref. 1) . 

Such conditions exist or are approached in an induction- stirred, 
melt- in- vacuum or liquid drop" in- space . 

This memorandum deals with the normal evaporation of binary 
alloys. In particular, we study the evaporative segregation pat- 
terns, i.e., the type and degree of enrichment or depletion of 
solute in the evaporating source, at different evaporation tempera 
tures and times. 
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DIFFERENTIAL EQUATION OF NORMAL EVAPORATION 


The appendix of Ref. 2 gives the exact solution for the normal 
evaporation of binary alloys. This equation relates the concentra- 
tion of the solute, m, at a time, t, when the mole fraction of 
the alloy remaining is F, as follows: 
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m / VI 
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U-V 


(1) 


where N and N q are, respectively, the number of moles of both 
solvent and solute at evaporating times t = t and t - 0; m is 
the initial molar concentration of the alloy; and U and V are, 
respectively, the evaporation rates of pure solute and solvent. 

For pure solute and solvent, respectively, these evaporation 
2 

rates in mol/cm /sec are 

A “B /T i 

U = K 10 U U (M u T) 2 (2) 

and 


ri D / x _ i_ 

V * K 10 (M y T) 2 (3) 

where K - 5.833 x 10 a, a = 1 for most metals (Ref. 1); M u 
and M v are, respectively, the molecular weights of the solute 
and solvent atoms; T is the evaporating temperature in degrees 
Kelvin; and A^ and or A and B v are evaporating con- 

stants .for the solute or solvent given, for example, by Ref. 1. 
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Differentiating Eq. (1) yields 


but 


dm 


(U - V)m(l - m) 

[mU + (1 - m)V]N aN 


dN - - [mU 4- (1 - m) V]A • dt 


( 4 ) 


(5) 


where A is the evaporating area of the alloy, assumed constant 
here. 

Substituting Eq. (5) into Eq. (4) results in 


dt = Gm a (l - m)^dm 


( 6 ) 


where 


G * 


N (1 - m ) 
_Q_ SI 


a+2 


A(U - V)m 


a-fl 


(7) 


a = (2V - U)/(U - V) (8) 

P = (V - 2U)/(U - V) (9) 
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SPECIAL CASE SOLUTIONS OF THE 


NORMAL EVAPORATION EQUATION 

Equations (6) -(9) allow us to determine the evaporating time, 
t, for an alloy to reach a specific solute concentration, m. Un- 
fortunately, these equations are not exactly solvable in the gen- 
eral case. All of the following special and important cases, (ex- 
cept Case V), however, are solvable in closed forms. 

Case I: The solute is much more evaporative than the solvent, 

i.e., U » V; or a — -1 and £ « -2. In this case 

G. dm 

i (10) 
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( 11 ) 


When m ^ m Q , the second term in the bracket is nearly 0. 

The evaporation time for this first case, t^, is then a logarith- 
mic function of m, as has been experimentally observed (see e.g.. 
Ref. 3) . 

Also, when m ~ m Q ~ 0, 

N m 

t-, ~ r® in - a (12) 

1 AU m v ' 
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Case IX: The solvent is much more evaporative than the 
solute, i.e., V » U; or a = -2 and p - -1. In this case 


dt 


G 2 dm 


^ (1 - m)m^ 


(13) 


where 


and 


m N 

r = — 2—Q. 
G 2 AV 


c 2 “ G 2 


(1 - m )m m - m 
in " £L ~r + tt 


m (1 - m) 
o x ' 


m m 
o 


(14) 


For dilute solutions (i.e., m ~ m i 0)> t 2 also becomes a 
logarithmic function of m, as has been observed. 

Case III: a = (2V - U)/(U - V) ■ 0, i.e., U = 2V and p = -3 


to — G. 


f 


-2 


-21 




(15) 


where 


G 3 


V 1 ~ m o> 

2AVm 


Case IV: P = (V - 2U)/(U - V) = 0, i.e., V « 2U and a = -3. 

In this case 



(16) 
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where 


G 4 " 


N m 2 
- o. o 


2AU(1 - m Q ) 


Equations (15) and (16) show that under some conditions, the 

evaporating time, t, is more adequately represented by linear 

-2 -2 

functions of (1 - m) or m , rather than by logarithmic func- 
tions of m. 

Case V: For relatively dilute alloys, i.e., m and Pm « 1, 

the following solution by series expansion can be obtained from 
Eq. (6): 


dt,. = Gm 
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(18) 


For computer calculations, it is desirable to know the ratio 
of the i*"* 1 term to the (i - 1) ^ term, thus 


T t P(a + l)(m 


a+i+1 a+i+l\ 
- m ) 

o / 


T ±-i 


i(a + i + 1)^ 


a+i a+i\ 


(19) 


m - m 




which is generally less than Pm/i or Pm Q /i. 
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Because pm « 1 and i constantly increases with each addi- 
tional term, this series converges rapidly unless p is very large, 
i.e., unless U = V, which leads to the following interesting case: 

Case VI; The solute and solvent are evaporating at equal rates, 
i.e., U = V. In this case, we would expect 

m = m Q for all t . (20) 

There is, then, no evaporative segregation, that is, there is 
neither solute enrichment nor depletion in the evaporating source. 

For any pair of solvent and solute, there is a unique tempera” 

ture, T , at which U = V and, hence, the alloy concentration 
3 

remains stable or constant. Equations (2) and (3) give: 

B - B 

T . U 3Z 

s \ - A v - 0.5 log(M u /M v ) ' 

Case VII: With extremely dilute alloys, i.e., m ~ m Q ~ 0, 

we have (Ref. 2) 


F « N/N q = (m/m o ) a+1 

Hence , 


( 22 ) 


t, = 


N r 

* _al 

AVl 


1 - (m/m 0 ) a+1 j 


(23) 
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EXAMPLES OF COMPUTATION 


• The Ni~Al System: 

As an example of the use of the various derived equations, 

the evaporation behavior of an alloy containing 8 percent by weight 

of A1 in Ni at the melting point of pure Ni (i.e., 1453°C) 

is computed. Here, the solute element (Al) is comparatively 

highly evaporative relative to the solvent (Ni) . Equation (11) , 

therefore, applies, and the time, tp to reach a final solute 

concentration m from a specified initial concentration, m Q , is 

directly proportional to N q A/U (in the G-^ constant) . Table 1 

gives the times to reach various final Al concentrations for one 

mole (53.66g) of the 8 percent Al in Ni alloy (m = 0.159) 

° 2 

evaporating at 1453 °C from its (supposedly constant) 10 cm 
surface. 


NORMAL EVAPORATION OF 8% BY WEIGHT ALUMINUM IN NICKEL AT 1453 °C 


Final Solute 

Concentration Alloy Remaining 


Evaporation 
Time (sec) 


Weight 

Fraction 

Mole 

Fraction 

Weight (g) 

Mole 


0.080 

0.1590 

53.66 

1.0000 

0.00 

0.079 

0.1573 

53.60 

0.9977 

14.08 

0.078 

0.1555 

53.53 

0.9955 

28.29 

0.077 

0.1536 

53.47 

0.9932 

42.64 

0.076 

0.1518 

53 .41 

0.9910 

57.12 

0.075 

0.1500 

53.34 

0.9887 

71.75 

0.074 

0.1481 

53.28 

0.9865 

86.54 

0.073 

0.1463 

53.22 

0.9842 

101.46 

0.072 

0 . 1444 

53.15 

0.9820 

116.55 

0.071 

0 . 1426 

53.09 

0.9798 

131.80 

0.070 

0.1407 

53.03 

0.9775 

147.22 
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• Iron Alloys: 

The evaporation behavior of binary iron alloys containing 20 
different solute elements has also been studied. Table 2, listing 

Table 2 

EQUI-EVAPORATIVE TEMPERATURES FOR 
TWENTY BINARY IRON ALLOYS 


Solute 

U/V at 
1600 °C 


T s , °C 

Solute 

U/V at 
1600 °C 

T , °C* 

S ■ * 

Cd 

3.765 x 

io 6 

13,540 

W 

4.412 

X 

io " 12 

-69,590 

Zn 

2.065 x 

10 6 

15,610 

c 

8.210 

X 

io " 8 

5,263 

Mg 

9.074 x 

10 5 
4 

19,610 

Mo 

1.879 

X 

io - 7 

-6 

-12,130 

Ca 

5.503 x 

10 

9,335 

Zr 

1.950 

X 

10 

-47,790 

Pb 

3.333 x 

10 3 

4,976 

B 

6.832 

X 

io ' 5 

9,510 

Mn 

1.070 x 

10 3 

20,740 

V 

3.804 

X 

io “ 3 

8,574 

Al 

45.65 


7,918 

Ti 

2.254 

X 

io ' 2 

34,650 

Cu 

12.25 


5,613 

Co 

4.362 

X 

io" 1 

4,318 

Sn 

9.978 


2,686 

Si 

5.238 

X 

io' 1 

2,825 

Cr 

3.158 


-215 

Ni 

5.897 

X 

io ' 1 

3,036 

the equi- 

-evaporative 

temperatures 

for these 

20 different 

alloy sys- 


terns, also gives the ratios of the solute evaporating rate, U, at 

1600 °C, to that of the solvent iron, V. At 1600 °C, 10 of these 

solute elements evaporate faster than the solvent (the three left 

columns) and 10 slower (the three right columns) . Moreover, these 

-12 

ratios vary widely over 18 decades, from 4.41 x 10 for the 
slowest evaporating, W, to 3.76 x 10 8 for the fastest evap- 
orating, Cd. Because of this wide variation in evaporating rates. 
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and because of the extreme sensitivities of the evaporating surface 
to unsuspected contaminants, predicted or experimental evaporating 
results cannot generally be very accurate. 

The equi-evaporative temperatures in iron alloys also vary 
widely. Binary iron alloys containing Cr, Zr, Mo, and W have 
no practical equi-evaporative temperatures. One can, therefore, 
always expect these alloys to change compositions continuously 
with the evaporating time. 

Table 3 shows the effect of evaporating temperature on the 
U/V ratios for four different solute elements Mg, Ca, Mn, and 
Al. In the range of 1500°C to 1900°C and beyond, increasing 

Table 3 


EFFECT 

OF SOLUTE 

ELEMENTS AND EVAPORATING TEMPERATURES 

ON 

THE VALUES 

OF U/V (x 1000) 

IN IRON 

ALLOYS 



Solutes 


Temp 

(°C) 

Mg 

Ca 

Mn 

Al 

1500 

2132.0 

118.2 

1.649 

0.06037 

1600 

907.4 

55.03 

1.070 

0.04565 

1700 

421.2 

27.68 

0.7261 

0.03552 

1800 

210.5 

14.88 

0.5110 

0.02831 

1900 

112.1 

8.466 

0.3719 

0.02304 


the evaporating temperatures always decreases the U/V ratios. 

This can also be seen from Table 2, as the equi-evaporative temper- 
atures shown for these four binary iron alloys are higher than 
1600 °C, at which temperature the four solute elements evaporate 
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much faster than the solvent iron. Table 3 also shows that for the 
same temperature variation, the more evaporative the solute element, 
the more percentage variation the U/V ratio. In the case of Mg, 
e.g., the U/V ratio decreases by about 20 times from 1500°C to 
1900 °C, whereas for Al, the same ratio decreases only by less 
than three times over the same temperature interval. 

Figure 1 shows the effect of solute elements and evaporating 
temperature, T, on the evaporating time t, for a given set of 

initial and final solute concentrations (i.e., m Q = 0.01 and 

m = 1 ppm) . For given m Q and m, the log t versus 1/T curves 
for these highly evaporative solute elements are approximately lin- 
ear and have positive slopes. This can be expected from Eq, (11) 
since, for given m Q and m, log t is linearly related to 

B/T - 0.5 log T, and since 0.5 log T is small relative to B/T 

within the evaporating temperature range studied. Thus, one can 
determine, from Fig. 1, the value of the elemental evaporating con- 
stant B, or heat of evaporation AH = 4.574 B (Ref. 1), for the 
solute elements by plotting log t versus 1/T and measuring the 
slope of the resultant, nearly straight lines. 

In Fig. 2, the log t versus 1/T relationships also appear 
nearly linear for all the four different initial concentrations 
(i.e., m Q = 10 \ 10 10 and 10 ^) of Al in Fe. This 

can also be seen from Eqs. (11) and (12). All these nearly straight 
lines have identical slopes, from which the heat of evaporation of 
pure Al can be evaluated. 

Figures 3 and 4 display the effect on the evaporating time, 
t, of final solute concentration, m, and either initial concen- 
tration m Q (in Fig. 3) or solute elements (in Fig. 4). 

Figure 5 indicates that Mg, Ca, and Mn in Fe alloys are 
so evaporative at 1600 °C that practically all of these elements 
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SEC 



1900 1800 1700 1600 1500 

EVAPT, °C 


SOLVENT: Fe 

m =0.01 
o 

m = 1 ppm 


Fig. 1 Evaporating Times for Different Dilute Iron Alloys 

to Change Concentration from = 0.01 to m = 1 ppm 
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EVAPt, 

SEC 



i i i i i 

1900 1700 1500 

EVAP T, °C 


Fig. 2 Evaporating Times for Dilute Iron Alloys Containing 
A1 and Starting at Four Different Initial 
Concentrations, to Reach a Final Concentration of 
m = 1 ppm 
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1 ( 

FINAL CONC, m 


Fig. 4 Effect of Solute Elements and Final Concentration 
on the Evaporating Times in Dilute Iron Alloys 
at 1600 °C 



1.0 r 


LOSS OF 
MATERIAL, 
-AN 


0.1 


0.01 



Mg, Ca, Mn 


0.001 


• SOLVENT = Fe 

• m =0.01 

o 

• T = 1600 °C 


0. 0001 1 , , , i i i — 

io " 2 io ’ 4 io ’ 6 io ' 8 


FINAL CONC, m 


Fig. 5 Evaporation Loss of Dilute Iron Alloys at 1600 °C 
As a Function of Solute Elements and Final Solute 
Concentrations 
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are removed by evaporation, without much evaporative loss of the 

solvent Fe atoms. On the other hand, A1 is comparatively less 

evaporative so that much of the solvent Fe atoms are evaporated 

off together with Al* To achieve a purification factor of 10 ^ 

-9 

(i.e., to m = 10 ) from m Q = 0.01, for example, the initial 

evaporating alloy must lose over 30 percent of its material. 
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ACCURACY OF PREDICTED RESULTS 


To check the validity of our derived equations for normal 
evaporation of binary alloys, the current literature has been 
searched. Several sets of alloy evaporation data have been found 
that are amenable to normal evaporation analysis. These sets in- 
clude the following: 

• The Fe-Ni and Ni"Cr Systems: 

An analysis has been made of the data by Obradovic et al. 

(Ref. 4) for Fe“Ni and Ni“Cr alloys evaporated at 1600 °C 

for different times under various ambient pressures . In the 

80% Ni-20% Cr case, evaporation started with N q = 0.7582 moles 

of the alloy having an evaporating area A = 10.0 cm^. The solute 

and solvent evaporating rates are, respectively, U_ = 3.479 x 10 

“6 2 

and = 6.386 x 10 mol/cm /sec. This is a case where the 

solute, Cr, is comparatively highly evaporative relative to the 
solvent, Ni. Hence, the evaporating time t-^ given in Eq. (11) 
applies, i.e.. 


t 



(1 - m Q )m m - m 

^ “ q ^ 1 " m ) + (1 " m Q ) (1 " m ) j 


where G-, = ~N (1 - m )/AU, and m and m are the initial and 
l o o o 

final molar solute concentrations in the alloy. 

Least square fits of the Obradovic data give nearly constant 

4 

values of observed G-, : ”1.785, “5.064, and -5.064 (x 10 sec) 

1 -3 

for ambient pressures of 1, 100, and 500 (x 10 torr) . The cal" 
culated values of G^, though also nearly constant, are, however, 
one order of magnitude smaller, indicating surface solute depletion. 
Similar analysis for the 60% Fe-40% Ni alloy evaporated with 


18 



9 — fi 

N = 0.7478 moles and A = 10.0 cm gives ll T . = 6.386 x 10 

° -5 2 Nl 

and Vp e = 1.081 x 10 mol/cm /sec. Here, the solvent Fe is 

much more evaporative than the solute Ni, and Eq. (14) applies. 

The values of observed evaporating coefficients G 2 = m Q N o /AV are 

2.469, 1.369, and 8.386 (x 10 sec), for ambient pressures of 

-3 

1, 100, and 500 (x 10 torr) , respectively. The calculated 
G 2 *s are again nearly constant but also an order of magnitude 
smaller than the observed G 2 *s, again indicating solvent deple- 
tion (or solute enrichment) at the surface. Details of these 
analyses are given in Ref. 5. 

• Beryllium Purification 

The kinetics of normal evaporation for beryllium have been 
quantitatively checked with actual results of beryllium purifica- 
tion during vacuum induction melting. Details are given in Ref. 6. 

In these tests, beryllium was crucible- free, induct ion -melted under 

*”6 

an ambient pressure of 10 torr. The actual temperature of the 
melt was not known, and the exact evaporating rates for the beryl- 
lium solvent and various solute elements cannot be computed. How- 
ever, because all solute elements (i.e., Fe, Cr, Mn, Ni, Si, Al, 
Mg, Cu, Zn, and Na) and the solvent. Be, were evaporating 
from the same or common liquid-gas interface of a fixed area for 
the same length of time, we can compute the values of P, defined 
as the product of th% evaporating, tiffife* t* and solvent evap- 
orating rate,, ¥ (i.m ¥ m ?§£«* § 

(for U » v) ; thus: £ » t¥ or tU, h §1 ¥§§ &g 

concentrations and final e©ne©n£rati©p§ (ip ppa) ©f ¥§£i8H§ 
solutes, together with tht value ©f §§|y§l §urf§g@ cgpgeptrp- 

tions, ratio of actual aurfsee t© hulk §@M§hi?§ki8fl§ 1 feffge- 
tive times to reach the final §©n§@n|r§|i©p§ $|g 

that the evaporating temperature wa§ 1110 8 §* l@ii8Wih§ S8H" 

elusions can be drawn from Tafel© k\ 



1 . 


The solute elements can be divided into three groups: 

Fe, Cr, Ni, Si, and Cu evaporate much more slowly 
than Be; Mg, Zn, and Na evaporate much more rapidly 
than Be; A1 evaporates at about the same rate as Be. 

2. The computed times to reach the final concentrations 
under normal evaporation conditions (complete liquid 
mixing) are fairly constant for the solute elements 
that evaporate much slower than Be (i.e . , Fe, Cr, 

4 

Ni, Si, and Cu) , being about 4 x 10 . 

3. The computed times to reach the final concentrations 
for the highly evaporative Mg, Zn, and Na are also 
fairly constant, but about four orders of magnitude 
smaller. 

4. After correction for limited liquid mixing (Ref. 7), 
the effective times to reach the final concentrations 
are much more constant, even between the groups of 
solute elements. In particular, the highly evapora- 
tive elements Mg, Zn, and Na have their results 
improved by several orders of magnitude. 

5. Limited liquid diffusion must, therefore, be con- 
sidered, particularly for the highly evaporative 
elements . 

6. The surface solute concentrations can be changed by 
up to six orders of magnitude, so that the effective 
evaporating rates can be similarly changed. 
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PURIFICATION KINETICS OF BERYLLIUM DURING CRUCIBLE -FREE VACUCM INDUCTION MELTING 
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DISCUSSION AND CONCLUSION 


These and other similar analyses suggest that the normal evap- 
oration equations presented in this Memorandum work surprisingly 
well for the several alloy systems studied so far. These equations 
generally give results accurately to within an order of magnitude 
of the correct values, perhaps even for highly concentrated solu- 
tions (e.g., 207 o Cr in Ni or 40% Ni in Fe) . Such accuracies 

are usually sufficient for many evaporation studies. 

In other more critical studies, however, refined evaporation 
analyses may be needed. In particular, the effect of limited dif- 
fusion and, hence, surface depletion or enrichment of the solute, 
must often be accounted for. This is true especially for the 
highly evaporative solutes such as, e.g.. Mg* Zn, and Na in 
the solvent Be. We already have results that confirm this con- 
clusion. Details of these results will be published in Ref. 7 and 
elsewhere . 
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Captions 


Figure 

1 Evaporating times for different dilute iron alloys to 
change concentration from m Q = 0.01 to m = 1 ppm 

2 Evaporating times for dilute iron alloys containing 

Al and starting at four different initial concentrations, 
to reach a final concentration of m = 1 ppm 

3 Effect of initial and final solute concentrations on the 
evaporating times in dilute iron alloys containing Al 
at 1600 °C 

4 Effect of solute elements and final concentration on 
the evaporating times in dilute iron alloys at 1600 °C 

Evaporation loss of dilute iron alloys at 1600 °C as a 
function of solute elements and final solute concentra- 
tions 
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